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Abstract 

In any scenario of the electroweak baryogenesis, the profile of the CP violating bubble 
wall, created at the first-order phase transition, plays an essential role. We attempt to 
determine it by solving the equations of motion for the scalars in the two-Higgs-doublet 
model at the transition temperature. According to the parameters in the potential, we 
found three solutions. Two of them smoothly connect the CP-violating broken phase and 
the symmetric phase, while the other connects CP-conserving vacua but violates CP in 
the intermediate region. We also estimate the chiral charge flux, which will be turned 
into the baryon density in the symmetric phase by the sphaleron process. 
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1 Introduction 



Since the fascinating proposal of the electroweak baryogenesis its various aspects have 
been investigated by many authors. For the baryogenesis to occur, the following conditions 
must be met: (1) The electroweak phase transition (EWPT) must be first order to realize 
a state out of equilibrium. (2) There had to be CP violation in the era of the EWPT. 
Besides these, to keep the generated baryons, (3) the sphaleron processes must decouple 
just after the EWPT. These impose some restrictions on the models of the electroweak 
theory. The conditions (1) and (3) give an upper bound on the lightest neutral Higgs 
particle, which is inconsistent with the present lower bound of the Higgs scalar in the 
minimal standard model. Further to have efficient CP violation, an extension of the 
Higgs sector would be needed. On this ground, the two-Higgs-doublet model, including 
the minimal supersymmetric standard model (MSSM), has been studied to estimate the 
generated baryon number. 

It is essential in the scenarios of the baryogenesis, spontaneous or charge transport, to 
know the profile of the CP violation near the expanding bubble wall created at the EWPT. 
In the literatures, however, some functional forms of the CP violation were assumed 
without any reasoning. They must be determined by the dynamics of the scalar fields near 
the EWPT. One may expect that, in the lowest order of the approximation, spacetime- 
varying CP violation would be governed by the classical equations of motion of the gauge- 
Higgs system, in which the Higgs potential is replaced with the effective potential at the 
transition temperature. This amounts to find the critical bubble, which would be a good 
approximation to an expanding bubble if the EWPT proceeds calmly. 

In this paper, we shall follow this line to obtain the functional form of the CP- violating 
phase in the two-Higgs-doublet model. In section 2, we derive the equation for the CP- 
violating phase assuming that the moduli of the Higgs scalars take the kink shape with 
the same width. In section 3, we enumerate possible boundary conditions of the phase for 
various choices of parameters in the potential. Next we present some numerical solutions 
and the chiral charge flux in section 4. The final section is devoted to discussions. 

2 The Equation for the Phase 

2.1 The Equations of Motion 

The system we concern is governed by the lagrangian, 

C = --F; v F a ^- l -B^+ J2 {D^if D^&i — V e ff(<&i, <J> 2 ; T), (2.1) 

4 4 t=l,2 
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where 



= (d, - ig T -Al{x) - i g -B,{x))^{x) 



We adopted V e ff at T near the transition temperature as the Higgs potential. We are 
now interested in the classical solution of the bubble-wall shape which mediates between 
the broken and symmetric phases. If the phase transition proceeds calmly, it will be valid 
to expect that the bubble wall grows keeping the profile of the critical bubble, which is 
determined by the static equations of motion. Further, when the bubble is spherically 
symmetric or is sufficiently macroscopic so that it is regarded as a planar object, the 
system is reduced to an effective one-dimensional one. In general, in 1 + 1-dimensional 
gauge theories, gauge fields have no dynamical degrees of freedom, that is, they are pure 
gauge. Here we assume that the gauge fields are written in the pure-gauge form: 

ig^A^x) = d,U 2 (x)U^(x), i 9 -B,(x) = dJJ^U^x), 

where U 2 and U\ are elements of SU(2) L and U (l)y, respectively. Since we can completely 
gauge away these gauge fields, we only need to consider the equations of motion for the 
Higgs fields. Assuming that U(l) em is not broken anywhere, the vacuum expectation 
values (VEVs) of the Higgs fields have the following form: 



<**(*)>= [ ±p t (x)e W > {x) < < = 1 ' 2 )- (2 ' 2) 



Now the equations of motion are 



8V tf 

-d'p^+p^d^d^x)-^- = 0, 

d»(p*(x)d»e i (x)) + ^ = o. 

Along with these equations, we have the "sourcelessness condition" , which arises from the 
requirement for the gauge fields to be pure-gauge type: 

f%(x)d li e 1 (x)+ f %(x)d ll 6 2 (x) = 0. 

Regarding the bubble wall as a static planar object, these equations are reduced to (taking 
z as the coordinate perpendicular to the wall) 

,d6i{z) . dOoiz) 
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For later convenience, let us change the variable z to y defined by 

1 



y 



(1 — tanh(az)) 



(2.6) 



where a 1 characterizes the width of the wall. In terms of this new coordinate, the above 
equations are written as 



Aa 2 y(l-y)-^- 
dy 



yiX-y) 



dp % {y) 

dy 



Aa 2 y(l-y)-^- 
dy 



yi 1 - y)p 2 i(y) 



, d y J 

d9 i (y) 



dV f 



eff 



dy 



Piw— i — + f4w 



dpi 

dVeff 

09 i 
d8 2 (y) 



0,(2.7) 
0,(2.8) 
0.(2.9) 



dy dy 

In order to solve these equations, one must know the explicit form of V e ff. Because of 
the gauge invariance, V e ff is a function of 81 — 82- From this fact, ( |2.8|) with % = 2 is 
automatically satisfied as long as pi and 8i satisfy ( |2.8| ) with i — 1 and (|~9|). 



2.2 Ansatz for the effective potential 

In general, it is difficult to solve the coupled equations ( [2.7|) and (|2.8j ) with the constraint 
( ^9|) . When the EWPT proceeds calmly, we expect that the modulus of the Higgs, Pi(z), 
has the shape of a kink, and that p\ and pi have the same order of width.[] Assuming 
that pi and P2 are the kink type of the same width but with different amplitudes, the 
problem is now to solve ( |2.8|) in the background of This assumption, in turn, restricts 
the form of the effective potential. We shall solve the equations for 9^ without specifying 
any model, but with the potential matching this assumption. 

We require that ( |2.7| ) has the kink-type solutions in the absence of CP violation (9 = 
9 1 - 2 = or tt); 

Pi{y) = Vi(l-y), 



(2.10) 



where 



Then 



vi = v cos P, v 2 = v sin j3. 



is 



4a 2 v iy (l-y)(l-2y) + 



dV ef f(pi,p 2 ,8 = or tt)) 



dpi 



0. 



(2.11) 



Pi=Vi(l-y) 



For this equation to be satisfied, the form of V e ff is somewhat restricted. Now let us 
determine V e ff in terms of a polynomial of p\ and p 2 - 

1 We implicitly assumed that both the VEVs acquire nonzero values at about the same temperature. 
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The most general tree-level Higgs potential is given by 

V = mi${$i + m2$2$2 + (ml$\$ 2 + h.c.) 

+ |A 1 (^t<^ 1 )2 + |a 2 (<^^ 2 ) 2 + A 3 (<^l^i)(«>^2) - A 4 (<^ + 1 ^ 2 )(^^i) (2.12) 

- jiA 5 ($I$ 2 ) 2 + [A 6 ($l$i) + A 7 ($t$ 2 )]($t$ 2 ) + h.cj , 

where m\ , m 2 , Xi, A 2 , A3, A4 G R and m|,A 5 ,A 6 ,A 7 G C, three of their phases are inde- 
pendent and yield the explicit CP violation. When all these parameters are real (i.e., no 
explicit CP violation), we have 

V o (pi,p 2 ,0) = ^mjpl + ^mlp 2 2 + ml Pl p 2 cos 9 + y Pi + yP2 ( 2 -13) 

A3 — A 4 2 2 A 5 2 2 (on\ ^63 a A 7 3 

+ 1 P1P2 - "J P1P2 cos(2#) - y P1P2 cos - y P1P2 cos 6. 

In the following we shall examine the spontaneously CP-violating and CP-conserving 
cases. The CP-conserving case is realized if 9 = or tt. It is sufficient to consider the 
former, since the latter is obtained by redefining one of the scalars. Without any CP 
violation, we have, in terms of pi, 

T/ 1 2 2 , 1 2 2 , 2 . Al 4 A2 4 

Vo = 2 m i p i + 2 m2P2 + m s^i^ 2 + y/^i + Y p2 

+ yPiP2-^(A6P? + A 7 p 2 )p lP 2, (2.14) 

where A3 = A3 — A4 — A5. 

In order to have kink solutions for P i, we need p 3 -terms with negative coefficients, 
which are expected to arise at finite temperature. Hence we adopt the following ansatz 
for the effective potential: 

V ef f(j> 1 ,p 2 ,e = 0) = V ( Pl ,p2,0) - (Ap\ + Bp\p 2 + Cp lP j + Dpi) . (2.15) 

Since we expect this to represent the effective potential with first-order phase transition, 
the origin (px, p 2 ) = (0, 0) and the point (p 1 , p 2 ) = (v cos fi, v sin f3) must be local minima. 
This condition amounts to 

det { tttt) > °« ^ > OT ^ > °» ( 2 - 16 ) 

\dpidpj) dpi dpi 

at each point. At (0, 0), this reduces to 

m\m\ — m\ > 0, and m\ > or m 2 > 0. (2-17) 

2 Although we do not use the high-temperature expansion, these terms may be considered to represent 
the effect which causes the first-order phase transition. 
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Requiring that ( J2.ll ) is satisfied with this potential, some of the parameters are expressed 
by the others. Now we assume cos (3 sin (3 7^ 0. This condition would be necessary to give 
masses to up- and down-type quarks when they couple to different Higgs doublets, as in 
the MSSM. Then the effective potential in the absence of CP violation is 



Veffipt, p 2 , 0) 



= (2a 2 - 



^3 tan (3)p\ + (2a 2 - ^m 2 cot (3)p\ + m\p x p 2 



-2 A cot f3 + D tan 2 (3 + 



4a 2 



;(3- 



1 



v sin (3 cos 2 (3 ' 



p\p2 



+ 



A cot 2 (3 - 2D tan (3 + 



4a 2 



;(3 



v cos (3 sin (3 ' 



P1P2 + Dpi 



A 2 



A, 



Al 4 A 2 4 A3 2 2 

+ y Pi + yP 2 + T PiP 2 



(2.18) 



+ 



3 A 8a 2 

-Ai cot (3 - — tan 3 /3 + A 3 tan (3 , 

2 2 v z sm p cos p 

Ai q . 3 , _ ~ _ 8a 2 

-— cot 3 /3+ -A 2 tan/3 + A 3 cot/3- — — 

2 2 v 2 sm p cos p 



(4- 



(4-- 



cos 2 (3 ' 
1 , 



sin (3' 



P\P2 



plpl 



and the condition (|2.17|) is written as 



4a 2 > 



;2.i9) 



sin (3 cos j3 

Note that although we use the same notations for the parameters in V e ff as those in Vq, 
their meanings are different. Those in V e ff contain radiative as well as finite-temperature 
corrections near Tq- 



Now let us introduce the CP-violating phase into (|2.18|) in a gauge-invariant manner. 
Comparing (|2.13|) and Q2.14j ) suggests that the phase 9 is introduced in V e ff of (|2.18|) as 



P1P2 

^p\p\ 

p\pi 
P\P% 



P1P2 cos#, 
A 5 p 2 p2 cos(26>), 
p\p 2 cos 8, 
Pxp\ cos 8. 



On the other hand, we have no principle to determine ^-dependence of p 3 -terms. Here we 
investigate the two possibilities: 

• no ^-dependence in the p 3 -terms, 

• p\p2 -> p\p2 cos 8, pipl p x pl cos 8. 
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Hence our ansatz for the effective potential is 



V eff {pi,p 2 ,6) 

(2a 2 ml tan f3)p\ + (2a 2 m 2 cot P)p\ + m\pxp 2 cos 9 

2 2 



-2A cot /3 + D tan 2 /? + 



4a 2 



r(3 



f sin /3 cos 2 p 1 ' 



p 2 p 2 (cos6>) 



+ 



Ai 



A cot 2 /3 - 2D tan /3 + 



4a 2 



(3-- 



A 



2 „4 



A, -A 



f cos /3 sin P 
A 



4 2 2 A 5 2 2 fnn\ 

P\P2 ~ ~rP\P2 cos(2#) 



4 ™ 4 
Ai cot (3 — tan 3 (3 + A 3 tan /3 — 



3, . „ A 2 



p x p\ (cos 9) + Dp 3 2 



(4- 



(2.20) 



v 2 sin /3 cos /? cos 2 /3 



+ 



A 3 

-^cot 3 /3 + -A 2 tan/3 + A 3 cot/3- - . ., 

2 2 it sin p cos p sin p 



P\Pl 

P1P2 \ cos 9. 



Here cos # in the p 3 -terms is unity in the case of the first possibility. 



2.3 Equations for the phases 

Once ^-dependence of V e ff is determined, one can derive the equations for 9{ from ( |2.8| 
In our case, the sourcelessness condition (12.91) is reduced to 



(2.21) 



cos 2 p— - + sin 2 p— - = 0. 
ay ay 

As noted above, we only need to solve the equation for 9%, as long as 9\ and 9 2 satisfy the 
sourcelessness condition. Since, from this condition, 9 2 (y) is written as 



2 (y) = -9t{y) cot 2 p + const. 



we have 



sin P 



(9i(y) + const. 



Noting that the derivative terms in the equation for 9\ are invariant under the shift of 9±, 
the constant in the r.h.s. of the above equation can be ignored. Putting 



o(y) 



sin P 



0i(y), 



(2.22) 



we have 



l - y) + y{l - y)(l - %)— — 



*/ 2 
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1 



m 2 sin (3 cos/3 - (A cos 3 [3 + D sin 3 /3 1 v(l - y) 



4a 2 sin (3 cos 2 /? 

?^ 2 ~ 16o 2 
+— (Ai cos 4 13 + A 2 sin 4 /3 + 2A 3 sin 2 f3 cos 2 /3 — )(1 - y) 



sin 6{y) 



+ ^-(1 - 2/) 2 sin%) cos%). 



(2.23) 



For later convenience, we denote the coefficients in the above equation as 



b = 



c = 



m% 



4a 2 sin f3 cos f3 ' 



32a 2 



(Ai cot 2 (3 + A 2 tan 2 (3 + 2A 3 ) - 



2 sin 2 (3 cos 2 (3 



d = 



8a 2 
4a 2 



( A 6 cot f3 + A 7 tan (3) , 



(2.24) 



4a 2 sin /3 cos 2 /3 



4a 2 

A cos 3 (3 + D sin 3 /3 



u / £ C 
+ 



4a 2 y sin (3 cos /3 / 

If there is no ^-dependence in the p 3 -terms in the V e ff, e = 0. Because of the condition 

(Hi), 

6> -1. 



(2.25) 



On the other hand, at (pi,p2) = (v cos f3, v sin /3) 



<9 2 K 



«9p 2 



kt 2 e - 12a 2 c + A 3 v 2 ) sin 2 f3, 



d 2 V e 



dp 2 
d 2 V e 



2 

eff 



dp 1 dp 2 



m 2 , — (8a 2 e — 12a 2 c + A 3 t> 2 ) cos 2 (3, 
m 2 + (8a 2 e — 12a 2 c + A 3 t> 2 ) sin (3 cos (3, 



so that 



\dpidpjj 

Thus the condition Q2.16|) is satisfied if 



16a 4 1 + b - 2e + 3c - 



hv 2 
4a 2 



6-2e + 3c> -1 + 



4a 2 



(2.26) 



Now ( |2.23| ) is written as 

y 2 (i-y) 2 ^+y(i-y)(i-^y) 



dd{y) 



d 



[b+c(l-y) 2 -e(l-y)}sm6(y) + -(l-y) 2 sm(26(y)). 

(2.27) 
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This is the equation that we shall examine in detail. One sees that 9(y) = mr with n£Z 
is the trivial solution. This equation is invariant under 9(y) i— > —6(y). This is because we 
have no explicit CP-violating terms in the potential. 

Before closing this section, we comment on the ansatz adopted. Although the kinks 
( p.lO| ) are solutions of (|2.11 ) with the potential ( P-18| ), they are no longer solutions of the 
coupled equations for pi and 9i with the potential ( |2.20|) . So that solutions to Q2.23Q will 
not be true solutions of the coupled equations. Further some of the parameters may be 
restricted to give a finite-energy solution for 9. We, however, expect that our solutions 
are not so different from the true solutions. This is because, as long as the p's have 
the kink shape, p's and 9's continuously rearrange themselves so as to take the minimal- 
energy configuration starting from our solutions. At the final stage, we should check that 
the physical quantities, such as the generated baryon number, are insensitive to small 
perturbations in the parameters of the potential. 



3 Asymptotic Behaviors of 6 



Before solving ( p.27|) , we investigate the asymptotic behaviors of the solutions. This will 
help us to find numerical solutions. Among possible solutions, we are concerned in those 
with finite energy density. When all the gauge fields are gauged away, the classical energy 
of the bubble wall is given only by the contribution of the scalars: 

E = /^x|E [^l(^)^)+V$t( x ).V$ i (x)]+y e// ($ 1 ,$ 2 ;T)J. 

For a static and planar bubble wall, the energy density per unit area is, in terms of pi 
and 9;, 



+ V eff (p 1 ,p 2 ,9) 



dy { ay(l -y) Yl 



i=l,2 



2ay(l - y) 



'dpj\ 2 / d9^ 

Av) P \dy , 

V eff (p u p 2 ,9) \ . 



(3.1) 



In the case of the kink-type wall treated in the previous section, for the energy density to 
be finite, we must have 



9'i(y) ~ y c "> with a > — 1 for y ~ 0, 
9[(y) ~ (1 - yf, with (3 > -2 for y ~ 1. 



(3.2) 
(3.3) 
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3.1 Asymptotic behavior in the broken phase 

The potential with the kink-type profile (|2.10 ) for pi is written as 



V eff (p u p 2 ,9) = 2aW(l-y) 2 

+4aV sin 2 /3cos 2 /3(l-?/) 2 {[o + c(l - yf - e(l - y)](l 

+±(l- y f(l- cos 2(9)1. 



COS I 



(3.4) 



Although in the symmetric phase the value of 9 is not determined because of V e ff(y = 
1) = 0, it will be allowed to take some specific values as shown later. On the other hand, 
in the broken phase (y — 0), 



v; 



eff\yo) 



4a V sin 2 f3 cos z (3 



(6 + c - e)(l - cos 9 ) + 7(1 - cos 26»o ) 



(3.5) 



where 9o = 9(0). If d = (A5 = 0), 9 = 2mr((2n + 1)tt) for positive (negative) b + c — e, 
i.e., no CP violation. When d 7^ 0, (|3.5| ) is written as 



V e ff(9o) = —2a v dsin /3cos (3 ^cos^o 
This implies that 



b + c — e 
d 



+ #o-indep. terms. 



(3.6) 



cos 9o 
9 



b + c — e 
-d ' 



for d < and \b + c — el < —a 1 , 



2n7r, for < —d <6 + c — e, ord>0 and 6 + c — e > 0, 
9 = (2n + l)7r, for 6 + c - e < d < 0, or d > and b + c - e < 0. 

Suppose that #(?/) is expanded as 

00 

0(y) = e o + f y Ea n y n (v>0,a ^0) 

n=0 



(3.7) 



at y ~ 0. The constraint z/ > matches the condition (|3.2|). When inserted in (|2.27|) , this 
yields 



IT ^ao + 



> + l) 2 ai -z/(2z/ + 3)a 



2/ 



n=2 

-(n + u- l)(2n + 2i/ + l)a n _i + (n + 1/ - 2)(n + v + l)a n _ 2 ] y n } 

OO "1 OO 

\2 



^ [6 + c(l-y) 2 -e(l- Z /)]cos^ 



OO 1 00 

.n=0 ^' n=0 



(3.f 
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H~(l-y) a cos(20 o ) 



2j2*ny n -^y 2u (J2 a ny n f + 

. n=0 ra=0 



+ [o + c(l-y) 2 -e(l-?/)] shift 
d 



1 oo 

z - n=0 



») 2 sin(20 ) 



r>2 OO 

1 - ^ 2v (E *ny n ) 2 + 

Z - n=0 



In order to have a nontrivial solution, 2v G Z. 

When #0 = n7r with n e Z, the lowest order terms give 

(6 + c — e) cos + d = v 2 . 



(3.9) 



The higher order terms of y will relate a n with n > 1 to the lower coefficients in a nonlinear 
way. 

When (9 7^ ri7r, v must be an integer for a nontrivial solution to exist. Then the lowest 
order terms yield 

o + c-e + dcosfl = 0. (3.10) 

For this to be satisfied, 

|6 + c-e|<|d|. (3.11) 

As we saw above, this 8q is energetically realized when \b + c — e| < — d with d < 0. For 
v = 1, 0(y)-terms of ( |3.8| ) lead to 

a o = [(0 + c — e) cos^o + dcos(26 ) )]ao — (2c — e) sin^o — dsin(2# ), 
which means, by use of ( |3.10| ), 

(2b - e) sin 6 



a 



1 + d sin 2 6*0 



(3.12) 



The higher order terms give relations among a n with n > 1 and the lower coefficients. 
For z/ > 2, 0(y)-terms of (|3.8| ) yield 26 — e = 0. When z/ = 2, 0(y 2 )-terms of (|3~8|) give 

(6 — e) sin #0 



o 



(3.13) 



4 + dsin 2 6> ' 

On the other hand, 0(y 2 )-terms for v > 3 lead to 6 = e, which means 6 = e = 0. We 
discard this case because of two reasons: One is because nonvanishing b oc m\ is needed to 
violate CP spontaneously when = as in the case of MSSM . The other is because m| 
in V e ff is induced in the presence of tree- level X^^j so that to have mj = is unnatural. 
Hence in this case, only v — 1 and v — 2 are allowed. Note that since ao = d'(0), the 
initial conditions for ( J2.27|j are completely fixed by the parameters of the potential. 
To summarize, the possible boundary values of 8q £1X6 clS follows: 
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(a) 8q = 2mr, if b + c — e > —d > 0, or d > and 6 + c — e > 0. The parameters must 
satisfy 6 + c — e + d = z/ 2 with 2z/ G Z. 

(b) # = (2n + l)ir, if b + c — e<d<0, or d>0 and 6 + c — e < 0. The parameters 
must satisfy —(b + c — e) + d = u 2 with 2v G Z. 

(c) cos 6*o = —(b + c — e)/d, if \b + c — e\ < —d and d < 0, for v = 1 or 2. 
Among these, we do not need to consider the case (b), as discussed in section 2. 

3.2 Asymptotic behavior in the symmetric phase 

Define ( = 1 — y and suppose that 9(C) is expanded as 

oo 

0( C) = e 1 + c »Y / b n ( n o>o, &o^o) (3.i4) 

?1=0 



at y ~ 1 (C ~ 0). In terms of ( |2.27|) is written as 



C 2 (l - + C(l - (3 - 4C)^P = (b + < 2 - eC) sin (9(C) + ^ sin(20(C)). (3.15) 

The condition 2/i 6 Z is required to have a nontrivial solution. The lowest order terms 
lead to 

&cos#i = Ma* + 2 ) ( 3 - 16 ) 

when #i = #(1) = ri7r, and 

6 = (3.17) 

and fi G Z when #x ^ nit. The higher order terms will give relations among the expansion 
coefficients. One may think that it is meaningless to ask the value of 9%, since in the 
symmetric phase the Higgs fields vanish so that CP is never violated in the Higgs sector. 
But what is important for the baryogenesis is how the phases of them tend to some value 
as the Higgs fields disappear. We shall not discuss the case of b = 0, i.e., m\ = for 
the reason stated above. Besides this, the first-order EWPT is realized when b > — 1 (see 
Q2.25| )). Hence the possible boundary values 9\ are; 

9 X = 2nir, b = //(// + 2) with 2fi G Z. (3.18) 
4 Numerical Analysis 

As we saw in the previous section, we are concerned with solutions satisfying the bound- 
ary conditions either (a) or (c), and (p. 18 ), according to the parameters in the effective 
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potential. We present numerical solutions satisfying each set of boundary conditions. One 
corresponds to the case in which CP is spontaneously violated in the broken phase, while 
the other has no CP violation in the broken phase. To find such solutions, we performed 
numerical analysis using the relaxation, as well as the shooting, algorithms, with the pa- 
rameters taken to satisfy the conditions ( 2.25 ) and (|2.26|) . We also evaluated the chiral 



charge flux with use of these numerical solutions, which is the basic quantity to generate 
the baryon number in the charge transport scenario 0. 

4.1 Solutions with spontaneous CP violation in the broken phase 

Although one may think that nonzero 8q is induced in some model at finite temperature H, 
we treat it as an input parameter. Once #0 is fixed, we can choose three of the four 
parameters (b, c, d, e) (see (c) in section 3). Because of ( |3.18|) , 6 = 5/4, 3, 21/4, 8.... 



As an example, we show the solution for the case of #0 = 0.002 and (6, c, e) = (3, 7, 7) in 
Fig. 1. This suggests that the real part of the VEV has the kink shape and the imaginary 
part can be regarded as a perturbation to it. Then the effects of the CP violation on 
the fermions scattered off the wall, which interact with the Higgs through the Yukawa 
coupling, can be treated by the perturbative method developed in |4[]. Recalling that 
Pi oc 1 — y for the kink- type profile, the behavior of 6(y) ~ 8q(1 — y) with small 6q 
in Fig. 1 implies that the imaginary part of the VEV is approximately proportional to 
the square of the real part, just as in the case studied in We shall not repeat the 
calculation of the chiral charge flux and note that A6 in [||] should be replaced with 8q 
here. 

Besides this solution, we also found a solution with 6> = 1 for (6, c, e) = (3,-1,0), 
shown in Fig. 2. Unlike the above solution, it no longer has a linear shape and the per- 
turbative method is not applicable. Such a solution would be realized when CP violation 
in the broken phase is enhanced by finite temperature effects and almost disappears at 
zero temperature as the case studied by Comelli, et a/.||. 

4.2 Solutions without CP violation in the broken phase 

In this case the parameters must satisfy (a). Then (p 1; p 2 , 6) = (v cos (3(1 — y), v sin j3(l — 
y),2nn) is an exact solution to the full equations of motion ( |2.7|) and Q2.8|) , which we 



refer to as the trivial solution. If there exists a nontrivial solution to (|2.27|) , it will open 
an interesting possibility that nonzero baryon asymmetry can be generated even in a 
CP-conserving theory in the broken phase. 



3 To O(6» ), 0{y) = (1 - yY with 2/i £ Z is an approximate solution to (|2.27|) for b = + 2), 
e = /i(2/i + 5) and d = + 3) — c < 0. 
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We found such a nontrivial solution for (b, c, d, e) = (3, 12.2, —2, 12.2), whose profile is 
plotted in Fig. 3. From this, the real and imaginary parts of the VEV are obtained and 
shown in Fig. 4, which suggest that the perturbative method is applicable. We calculated 
the chiral charge flux for various choices of the fermion mass m and wall width. The 
results are summarized as the contour plot in Fig. 5. This shows that the chiral charge flux 
is comparable to those studied in ||, so that it could generate sufficient baryon number 
for the thin wall case, taking into account the enhancement of forward scattering. 

The energy density per unit area of this solution is 

AS = £ - £\ e=0 = -2.056 x 1(T 3 av 2 sin 2 (3 cos 2 /3, (4.1) 

where £ is defined by ( |3.1| ) and £\ e=0 = av 2 /3. Since this solution is not a solution to the 
full equations of motion, the energy of the true solution may be lower than this value. We 
found that the solution is stable under perturbation of the parameters, so that we expect 
that this type of solutions will exist even if we do not impose the kink-type profile for p^. 

For such a nontrivial solution to exist, the CP- violating state must be favored in the 
intermediate region between the broken and symmetric phases. From ( |3.4j ), the effective 
potential along the kink is 

V eff ((,6) = 2aV{C 2 (l-C) 2 (4.2) 

-d( 4 sin 2 (3 cos 2 /3 



2 >\ 2 



' b + cC-eC Y , b + cC-eC 

cos 9 H — — 1 



d( 2 J V ^C 2 / 

where ( = 1 — y. Just as the boundary value 9q, 6(y) can take a value other than nir, 
when d < and there is a region in < ( < 1 such that |(6 + c( 2 — e()/(d( 2 )\ < 1. 
Further minimization with respect to ( in such a region will give the local minimum of 
V e ff. Since this new minimum appears on coupling the phase 9, the expected structure 
of V e ff(p\, p2, 6) is somehow modified compared with VJ=//(pi, p2, 0), which has degenerate 
local minima at (pi, p 2 , 0) = (0, 0, 0) and {y cos (3, v sin /3, 0) reflecting the first-order nature 
of the EWPT. Hence we must require that the new minimum with 9 ^ mr should not 
be so deep, otherwise it drastically changes the dynamics of the EWPT. Here we shall 
not investigate this condition further, but only assure that the new minimum is not lower 
than (0,0,0). For the above set of the parameters and (3 = 7r/4, for which the negative 
contribution is maximal, V e ff as a function of (y, 9) is shown in Fig. 6. The new minimum 
of V e ff is about 0.0069a 2 t> 2 and the height of the saddle point between it and the origin 
is about 0.038a 2 t> 2 , which is much lower than the maximum 0.125a 2 v 2 along 9 = 0. Since 
the barrier between the new minimum and the origin is not so high, we expect that the 
nature of the EWPT is not essentially altered. The contour plot of V e ff also shows that 
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the numerical solution for 9(y) goes around the maximum at (1/2,0). The value of V e ff 
where the maximum of 9(y) is reached is about 0.116a 2 u 2 , which is just below the CP- 
conserving maximum. This fact suggests that this type of solution may exist irrespective 
of the depth of the new local minimum. 

5 Discussions 

Based on a rather general effective potential in the two-Higgs-doublet model at the tran- 
sition temperature, which exhibits first-order EWPT, we have classified possible finite- 
energy bubble wall solutions with CP violation. We have found two types of numerical 
solutions, which are characterized by the CP- violating angle in the broken phase, 9q. 

One of them is the solution in the case where CP is spontaneously violated in the 
broken phase and will be the lowest-energy solution in that case. It smoothly mediates 
between the CP-violating vacuum in the broken phase and the CP-conserving one in the 
symmetric phase. For sufficiently small 9q, the profile of the bubble wall is similar to that 
used in the literatures to estimate generated baryon numbers, while for 9q = 0(1), 9(y) is 
no longer proportional to the kink so that the previous estimations of generated baryon 
number with large 9, by use of some presumed profile, should be revised. 

The other connects the CP-conserving vacua in both phases. When CP is not violated 
in both vacua, there always exists the trivial solution with 9 = all along. The new 
solution has nonvanishing 9 within the bubble wall and has lower energy than the trivial 
one, as shown in (|4.1| ). For the critical bubble of radius Rc, this CP- violating bubble will 
be nucleated with more probability than the trivial one by the factor 



According to the estimation in the massless two-Higgs-doublet model 0, the radius of 



the critical bubble is given by ^j3Fc/(4:7iav 2 ), where Fq is the free energy of the critical 

bubble and is found to be about 145T. Then the exponent in (|5.1|) is 0.89 sin 2 ft cos 2 ft, 

irrespective of the wall width. This solution will provide a new possibility to generate 

baryon number even within the framework of a CP-conserving theoryfj 

We estimated the generated baryon number with the use of each solution, based on 

the charge transport mechanism. With the enhancement of forward scattering, we could 

have sufficient baryon number for thin wall case. Although we did not work out on the 

spontaneous scenario, we expect that sufficient baryon number could be obtained if the 

effects of diffusion are considered 0. 

4 Strictly speaking, one needs a CP-odd term in the energy to have net baryon number by solving the 
degeneracy between the solutions with 8(y) and —9(y). 




(5.1) 
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We found that each type of the solutions are stable under perturbation of the parame- 
ters in V e ff, as long as they satisfy the required conditions. Hence these kinds of solutions 
will exist in wider class of models. It will be interesting to find a realistic model which 
admits such a solution. 

When we almost completed this work, we noticed a work done by J. M. Cline et ai, in 
which they found solutions with 9 ^ (|#o| <^ 1)H)- They incorporated the case of the 
explicit CP breaking, besides the spontaneous case, and their equation for 9 corresponds 
to that with c = e = in our analysis. 
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Fig.l: The numerical solution of 8(y) for 8(0) = 0.002 and 8(1) = 0. The parameters are 
(&,c,e) = (3,7, 7). 

Fig.2: The numerical solution of 8(y) for 8(0) = 1 and 8(1) = 0. The parameters are 
(b,c,e) = (3,-1,0). 

Fig.3: The numerical solution of 8(y) for 8(0) = 8(1) = 0. The parameters are (b, c, d, e) = 
(3,12.2,-2,12.2). 

Fig. 4: The profile of the bubble wall corresponding to the solution in Fig.3 with x = az. 
The solid line represents the kink, which is the absolute value of the VEV with the 
maximum being normalized to unity. The dashed line and dashed-dotted line are 
the real and imaginary parts of it, respectively. 

Fig.5: Contour plot of the chiral charge flux, normalized as Fq/ (uT 3 (Q l—Q r)) , for u = 0.1 
and T = lOOGeV, where u is the wall velocity and Fq is defined in 

Fig.6: Contour plot of V e ff/(a 2 v 2 ) along the kink, as a function of y and 8, for the param- 
eters used to find the solution in Fig.3 and for f3 = it /A. 
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